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We first give a short historical overview with some key facts of massive star population syn¬ 
thesis with binaries. We then discuss binary population codes and focus on two ingredients 
which are important for massive star population synthesis and which may be different in differ¬ 
ent codes. Population simulations with binaries is the third part where we consider the initial 
massive binary frequency, the RSG/WR and WC/WN and SNII/SNIbc number ratio’s, the 
probable initial rotational velocity distribution of massive stars. 


1 Introduction 

It is generally accepted that most of the Wolf-Rayet 
stars are massive hydrogen deficient core helium 
burning stars. Some WNL stars (using the origi¬ 
nal abbreviation-definition of Vanbeveren and Conti, 
1980) may be core hydrogen burning objects. They 
are not considered in the present paper. In the late 
sixties and seventies Roche lobe overflow (RLOF) 
in binaries was considered as a most plausible pro¬ 
cess capable to remove the hydrogen rich layers of 
a star. Furthermore, at the massive star conference 
in 1971 in Buenos Aires, Kuhi (1973) presented sta¬ 
tistical arguments to conclude that all galactic WR 
stars may be binary components. At the same time 
large observational data sets became available of X- 
ray binaries and together with the WR argument it 
made flourish massive close binary evolution. Some 
protagonists are B. Paczynski, E. van den Heuvel, 
I. Iben, A. Tutukov, L. Yungelson, C. De Loore and 
the Brussels gang. Interested readers may find many 
massive binary evolution studies published in this 
period. 

The early seventies are also characterised by im¬ 
portant breakthroughs in the study of stellar winds 
of massive stars (e.g., Castor et ah, 1975) and the 
question was raised whether or not WR stars could 
form via massive single star evolution, single stars 
that lose their hydrogen rich envelope by stellar 
winds (the Conti scenario, Conti, 1976). Chios! et 
al. (1978) were among the first to demonstrate that 
this is indeed possible, but it still had to be shown 
that WR single stars exist (remember Kuhi, 1971). 
Vanbeveren and Conti (1980) reconsidered the galac¬ 
tic census of WR binaries. They concluded that the 
Kuhi statistics is biased and that the real galactic 
WR-fOB frequency is no more than 30-40% a per¬ 
centage that still holds today. 

Since 1980 large evolutionary data sets of massive 
single stars and binaries were calculated and with 
these data sets in combination with a description of 
the different processes that govern binary evolution 
(e.g., the effect of the supernova explosion on or¬ 
bital parameters, the treatment of binary mergers, 
the treatment of the effect on orbital parameters of 
mass loss from the system during the RLOF, etc.) it 
became possible to predict theoretically how a mas¬ 


sive star population would look like. It was real¬ 
ized that a comparison with the observed population 
yields important information in order to understand 
the physics of massive star/binary evolution. 

Meurs and van den Heuvel (1989) are probably 
among the first authors to study in some detail the 
massive star population including binaries. The au¬ 
thors focussed on the number of evolved early type 
close binaries in the Calaxy, in particular the X-ray 
binary population. 

Portegies Zwart and Verbunt (1996) introduced 
the skeleton of what would later on become the pop¬ 
ulation code SeBa used in dynamical N-body simula¬ 
tions. In 1996 they used it to discuss the population 
of massive binaries with compact companions. 

The investigation of the effects of binaries on the 
population of 0-type and WR-type stars also started 
by the end of the previous millennium (Dalton and 
Sarazin, 1995; Vanbeveren, 1995; Vanbeveren et ah, 
1997; De Donder et ah, 1997). A detailed descrip¬ 
tion of the Brussels massive star population code 
was given in Vanbeveren et al. (1998a = Paper 
I) (see also Vanbeveren et al., 1998b for an ex¬ 
tended review) and it was applied in order to predict 
the galactic number ratios WR/0, WR-bOB/WR, 
WC/WN, O and WR stars with a compact compan¬ 
ion, O and WR runaway stars, O and WR single 
stars but with a binary origin, O and WR single 
stars originating from a merged binary, the number 
of real O and WR single stars, etc. Furthermore, De 
Donder and Vanbeveren (2004) studied the effects of 
binaries on the chemical evolution of galaxies and it 
was therefore necessary in order to extend the code 
so that it was capable to predict the temporal evo¬ 
lution of the massive star/binary population as a 
function of metallicity Z. 

Since 1998 other research groups wrote massive 
star-bbinary population codes with various degree 
of sophistication, e.g., Nelemans et al. (2001) and 
Toonen et al. (2012) substantially updated SeBa, 
Izzard et al. (2004, Binary_c) simulated the pop¬ 
ulation of core collapse supernovae and gamma-ray 
bursts, Belczynski et al. (2008, Startrack) focussed 
on relativistic binaries with an application to future 
gravitational wave detectors. Eldridge et al. (2008) 
also presented simulations of massive star popula¬ 
tions as a function of Z using an extensive grid of 
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single and binary evolutionary computations. 

Spectral synthesis is a powerful tool in order to 
investigate extragalactic young massive starburst re¬ 
gions in general, the O and WR spectral features 
in these starbursts in particular. Starburst99 (Lei- 
therer et ah, 1999) is extremely useful for galax¬ 
ies with active star formation but it has to be 
kept in mind that it does not account for the pres¬ 
ence/evolution of massive binaries. The effect of bi¬ 
naries on the spectral synthesis of young starbursts 
was a main research topic in Brussels (Van Bever et 
ah, 1999; Van Bever and Vanbeveren, 2000, 2003; 
Belkus et ah, 2003) and repeated by Eldridge and 
Stanway (2009) who essentially confirmed the Brus¬ 
sels results. 

De Mink et al. (2013, 2014) introduced an impor¬ 
tant parameter in massive star population studies: 
rotation. They showed that the observed distribu¬ 
tion of rotation velocities in 0-type stars can be ex¬ 
plained entirely by the process of RLOF and the 
possible spin-up of mass gainers/binary mergers. 

In the line of the work of De Mink et al. a study of 
the rotation velocities of the 0-type companions in 
WR-bO binaries may be most illuminating and may 
help to answer the question if RLOF and mass trans¬ 
fer happened in WR-l-0 progenitors. Mike Shara, 
Tony Moffat, Gregor Rauw, Dany Vanbeveren et al. 
just started an observational project aiming at deter¬ 
mining these rotational velocities in as many WR-l-0 
binaries as possible. We invite everybody interested 
in joining the et ah. 

2 Massive single star + binary 
population synthesis codes 

The Brussels population code has been described in 
detail in the list of papers given in the introduction 
and very recent updates are discussed in Mennekens 
and Vanbeveren (2014). Rather than repeating all 
the basic ingredients (once more) we prefer to high¬ 
light two ingredients which may be different in dif¬ 
ferent codes: the initial-final mass relation and the 
LBV/RSG stellar wind mass loss rate formalism. 

2.1 The initial-final mass relationship 

Some codes (e.g., Nelemans et ah, 2001; Izzard et ah, 
2004; Belczynski et ah, 2008; De Mink et al., 2013) 
use single star interpolation formulae proposed by 
Hurley et al. (2000) and binary-evolutionary algo¬ 
rithms described by Tout et al. (1997) and Hurley 
et al. (2002). Other codes (e.g., Eldridge et ah, 
2008; the Brussels code) use an extended library of 
detailed binary evolutionary computations. The lat¬ 
ter 2 codes and the codes based on the algorithmic 
method predict massive star populations that differ 
mainly in the absolute number of mass gainers of 


interacting-binaries but the differences are not crit¬ 
ical and do not affect overall conclusions made in 
the papers cited above. A much more severe differ¬ 
ence is related to the initial mass-final mass relation 
ship. As was discussed by Mennekens and Vanbev¬ 
eren (2014) for stars with initial mass > 20 Mq the 
final masses are significantly larger in Hurley et al. 
based codes than in the Brussels code, possibly due 
to different stellar wind mass loss rate formalisms 
and/or alternative convective core overshooting pre¬ 
scriptions during core hydrogen burning. Note that 
a similar effect is visible in the intermediate mass 
range (Toonen et ah, 2014). Unfortunately, this dif¬ 
ference plays a critical role for the predicted popula¬ 
tion of binaries with at least one compact compan¬ 
ion, even more for systems consisting of two compact 
stars, and as shown by Mennekens and Vanbeveren 
(2014) also for the predictions of the detection rates 
of gravitational wave observatories. 

2.2 The stellar wind mass loss rate 
formalism during the LBV and the 
RSG phase 

LBV. The Brussels code adopts the LBV scenario 
of massive binaries as has been introduced by Van¬ 
beveren (1991). It states that the LBV phase is 
a common evolutionary phase of the most massive 
stars and that the LBV mass loss rate suppresses 
the RLOF/common envelope phase in case Br/case 
Be/ case C binaries when the mass loser has a mass 
higher than «40 Mq (see also Mennekens and Van¬ 
beveren, 2014, for a recent argumentation). How¬ 
ever, when this mass loser is a member of a case A 
binary, the star will lose most of its hydrogen rich 
layers due to RLOF prior to the LBV phase and 
it is tempting then to assume that the LBV phase 
does not happen. A case A binary in this high mass 
range then follows a more canonical evolutionary sce¬ 
nario where the case A RLOF is followed by case Br 
RLOF. Mennekens and Vanbeveren (2014) demon¬ 
strated that the way how LBV mass loss is imple¬ 
mented in population codes has an enormous effect 
on the predicted merger rates of double compact star 
binaries (primarily double black hole binaries are af¬ 
fected) and thus also on the predicted detection rates 
of forthcoming advanced LIGO detectors. 

Note that the story would be completely different 
if the LBV phenomenon would appear to be related 
to massive binary mergers. 

RSG. Since the early days when scientist started to 
investigate the effects of stellar winds on massive 
star evolution, very conservative formalisms were 
used in order to study the effects of the RSG wind. 
Most common was the formalism proposed by de 
Jager et al. (1988). However, based on observations 
of Feast (1992), Bressan (1994) concluded that the 
mass loss rate during the RSG phase of an LMC 
20 M© star could be a factor 10 larger than pre- 
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dieted by the de Jager et al. formalism. Vanbeveren 
(1995) was among the hrst to investigate the effect 
of these larger rates on the evolution of 20 Mo - 25 
Mq single stars. It was concluded that as a con¬ 
sequence of RSG mass loss a 20-25 M© single star 
may become a WR star. A more throughout discus¬ 
sion of evolutionary computations of massive single 
stars with the alternative RSG rates and the effect 
on the overall WR population was presented in Pa¬ 
per I. Salasnich et al. (1999) also re-investigated the 
effect of new RSG mass loss rates on the evolution 
of massive stars and essentially arrived at similar 
conclusions. Note that these alternative RSG rates 
also significantly affect the evolution of case G bi¬ 
naries, the RSG scenario as it was described in the 
Vanbeveren et al. papers which states that RSG 
mass loss may suppress the RLOF in Case C mas¬ 
sive binaries. About 15-20 years after the papers of 
Bressan and Vanbeveren larger RSG rates were also 
implemented in the Geneva single star evolutionary 
code (Ekstrom et al., 2012; Meynet et al., 2015). 
In the latter paper it was concluded that enhanced 
mass-loss rates during the RSG phase have little im¬ 
pact on the WR population, contrary to the simula¬ 
tions made in Brussels. The difference between the 
Geneva and Brussels results is most probably due to 
the post-RSG mass loss rates used in both codes. To 
illustrate let us consider a 20 Mq star. When due to 
the larger RSG rates this star has lost about 10 Mq 
during the RSG phase, it leaves the RSG region and 
starts evolving to the blue part of the HR-diagram. 
The Geneva code calculates the further evolution by 
using blue supergiant mass loss rates (Vink-rates) 
and the star never becomes a WR star. However, 
although the star still has a rather high hydrogen 
content in its atmosphere at the moment it leaves 
the RSG phase (typically X^tm ~ 0.5), the models 
have an internal structure similar to WNL stars. In 
Brussels we therefore decided to continue the further 
post-RSG evolution by using typical WNL mass loss 
rates rather than blue supergiant rates. As a con¬ 
sequence the star loses its remaining hydrogen rich 
layers, becomes a WNE (also here we use the orig¬ 
inal nomenclature of Vanbeveren and Conti, 1980) 
and eventually a WC star. The discussion of post- 
RSG mass loss rates remains open but at least with 
the Brussels suggestion it is possible to explain the 
low luminosity WC stars as observed by Sander et 
al. (2012). In section 3.2 we will add additional 
support for higher RSG mass loss rates. 

3 Massive single star + binary 
population synthesis 
simulations 

Note that in this section we only consider population 
simulations where binaries are included. 


3.1 The initial massive binary frequency 

First, it is important to realize that accounting for 
all the physical processes that determine the evolu¬ 
tion of binaries, the massive binary frequency (in a 
population of stars where star formation has been 
going on for at least a few million years) is smaller 
that the binary frequency at birth (on the ZAMS). 
In all the population simulations that we published 
since Paper I it is assumed that the massive binary 
frequency (binaries with initial period < 10 yr) at 
birth f > 0.7. This latter value is based on the fol¬ 
lowing argumentation. By studying a sample of 67 
bright 0-type stars Garmany et al. (1980) concluded 
that 33% (±13%) are primary of a close binary with 
mass ratio > 0.2 and period P < 100 days. As dis¬ 
cussed in Paper I a population of 0-type stars in a 
field of continuous star formation consists of real sin¬ 
gle stars, interacting binaries with periods up to 10 
years, mergers looking like singles, post-supernova 
single 0-type mass gainers, etc.. Therefore, to re¬ 
cover the results of Garmany et al. (1980) with a 
population synthesis simulation we had to start with 
an initial binary frequency f > 0.7. 

As was outlined in the introduction, the observed 
WR±OB binary frequency (in the Galaxy) seems to 
be not larger than 30-40%. But also a population of 
WR stars in a region where star formation is con¬ 
tinuous consists of real WR single stars, WR±OB 
binaries, WR stars resulting from binary mergers, 
single WR stars resulting from post-SN single OB- 
type mass gainers etc. and also here we had to start 
with a very high initial binary frequency in order to 
explain the observed WR±OB frequency. 

Recent observations of 0-type stars in young clus¬ 
ters seem to confirm a high massive binary fraction 
(Sana et ah, 2008, 2009, 2011, 2012; Rauw et al., 
2009) and they propose a value > 50%. A similar 
exercise as the one made in Paper I (and summa¬ 
rized above) was done by De Mink et al. (2014) but 
using the observations discussed by H. Sana et al.. 
Also De Mink et al. concluded that in order to re¬ 
cover the observed > 50 % massive binary frequency 
of H. Sana and co-workers, one has to start with an 
initial binary frequency > 70%. 

3.2 the RSG/WR number ratio as 
function of Z 

RSG winds significantly affect the RSG-timescale of 
a massive star and depending on the post-RSG mass 
loss formalism (see section 2.2), they also signifi¬ 
cantly affect the WR-timescale. The RSG (and post- 
RSG) mass loss formalism therefore significantly af¬ 
fects the predicted RSG/WR number ratio. Remark 
that our population code is the only single star ± 
binary code that includes the effects of the alterna¬ 
tive RSG mass loss rates as discussed above. Fig. 1 
compares predicted and observed RSG/WR number 
ratios as function of Z. The observations come from 
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Massey (2003) with updates as reviewed by Massey 
et al. (2013). The predictions holds for a popula¬ 
tion with a binary frequency at birth = 70%. The 
dashed line is based on the simulation of Eldridge et 
al. (2008), the full lines are the Brussels predictions 
(predictions depend on parameters who’s values are 
subject to some uncertainty; varying the values of 
these parameters yields a maximum and a minimum 
RSG/WR number ratio, resp. the upper line and the 
lower line in Fig. 1). The main difference between 
the simulation of Eldridge et al. and the Brussels one 
is the RSG mass loss formalism and we are inclined 
to conclude that the prediction with the alternative 
(higher) RSG rates fits the observations better, sup¬ 
porting the conclusion of section 2.2. 



Fig. 1 : A comparison between the observed and pre¬ 
dicted RSG/WR number ratio as function of Z. The 
observations are those reviewed by Massey (2003) and 
the arrows indicate updates as discussed by Massey et 
al. (2013). The dashed line is the prediction proposed 
by Eldridge et al. (2008), the full lines are the max¬ 
imum/minimum (see text) predictions made with the 
Brussels code. 


3.3 The WC/WN number ratio as 
function of Z 

L. Smith (1973) argued that metallicity might be 
responsible for the relative absence of WGs in the 
Magellanic Clouds, but without understanding the 
physical mechanism. Vanbeveren and Conti (1980) 
were among the first to link the effect of Z on stel¬ 
lar wind mass loss and the WC/WN-Z dependence. 
Detailed population simulations (with binaries) were 
presented by Vanbeveren et al. (2007) and Eldridge 
et al. (2008) and compared to observations. It was 
concluded that correspondence is rather satisfactory. 
In both studies the observations were those discussed 
by Massey (2003). It can readily be checked that 
recent updates (Neugent et ah, 2013) do not signif¬ 
icantly change the overall conclusions made in the 
two population studies cited above. 


3.4 The SNII and SNIbc population 

Studies of the effect of binaries on the population 
of SNII and SNIbc are numerous and go back to 
the very beginning of massive binary evolution re¬ 
search. Interested readers may consider Tutukov 
et al. (1992), Podsiadlowski et al. (1992), Joss 
et al. (1992), De Bonder and Vanbeveren (1998, 
2003, 2004), Belczynski et al. (2002) and references 
therein. Some more recent work essentially confirms 
the earlier results. One of the conclusions is that 
most of the progenitors of SN Ibc are massive binary 
components with an initial mass > 10 Mq, e.g. most 
of the progenitors of SN Ibc do not have an initial 
mass > 25-30 Mq and thus most of the progenitors 
of SN Ibc are not WR stars. Moreover, since the 
population of SN II and SN Ibc depends so much on 
the massive interacting binary population, one may 
wonder whether SN-population differences between 
different types of galaxies may reflect differences in 
the population of these massive binaries. 

De Bonder and Vanbeveren (1998) compared the 
overall (cosmological) observed SN II and SN Ibc 
population with population simulations and it was 
concluded that the overall cosmological massive in¬ 
teracting binary frequency should be about 50%. 

3.5 The initial rotational velocity 
distribution 

The observed rotational velocities of 0-type stars in 
the Galaxy has been discussed by Gonti and Ebbets 
(1977), Penny (1996), and in Paper I. An analysis of 
these observations in terms of massive single star and 
close binary evolution was presented in Vanbeveren 
(2009). It was concluded that a majority of mas¬ 
sive 0-type stars are born as relatively slow rotators 
with an average < 200 km/s rather than the 300 
km/s used by the Geneva team. The rotation veloc¬ 
ity distribution proposed in the papers cited above 
shows that there is a significant group of rapidly- 
rotating 0-type stars but many of the stars in this 
latter group are runaways with a peculiar space ve¬ 
locity > 30 km/s. This means that many of them do 
not have a canonical single star history but are the 
product of binary evolution (binary mergers, spun- 
up binary mass gainers, mergers due to dynamical 
interaction in dense clusters). The latter paper then 
suggested the following: if one asks whether or not 
rotation is important for stellar evolution, the an¬ 
swer is yes but perhaps mainly in the framework of 
binary evolution (rapidly-rotating mass gainers and 
binary mergers) or in the framework of dynamics in 
dense clusters where stars collide, merge and become 
rapid rotators. 

Within the VLT-Flames Tarantula survey 
Ramirez-Agudelo et al. (2013) investigated the ro¬ 
tational velocities of the 0-type stars in 30 Dor and 
obtained a distribution which is very similar as the 
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one proposed for the Galaxy in the papers listed 
above. 

De Mink et al. (2013) implemented rotation and 
the evolution of rotation in a population code of sin¬ 
gle stars and of close binaries. They concluded that 
starting with an initial population of slowly rotating 
massive stars (average velocity = 100 km/s), the ob¬ 
served rotational velocity distribution (the one dis¬ 
cussed in the papers cited above) can be recovered by 
properly accounting for all processes that affect the 
rotation in single star and binary evolution. There¬ 
fore (in line with the suggestion of Vanbeveren, 2009 
cited above) it can not be excluded that most mas¬ 
sive 0-type stars are born as slow rotators, much 
slower than the average value adopted by the Geneva 
team in their standard evolutionary calculations. If 
this is true then one may be inclined to conclude that 
the overall evolution of massive single stars and of 
most of the binary components prior to the onset of 
RLOF hardly depends on rotation. 

4 Conclusions 

We like to end this paper with an advice and an 
overall conclusion. 

Advice (not only for young scientists): before start¬ 
ing a research topic try to get a literature overview 
that is as complete as possible and do not forget 
that also in the previous millenium interesting stud¬ 
ies have been published. 

Conclusion: a theoretical population study of mas¬ 
sive stars where binaries are ignored may have an 
academic value but may be far from reality. 
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